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Abstract — We consider cooperative relay communication in a 
fading channel environment under the Orthogonal Amplify and 
Forward (OAF) and Orthogonal and Non-Orthogonal Selection 
Decode and Forward (OSDF and NSDF) protocols. 

For all these protocols, we compute the Diversity-Multiplexing 
Gain Tradeoff (DMT). We construct DMT optimal codes for the 
protocols which are sphere decodable and, in certain cases, incur 
minimum possible delay. 

Our results establish that the DMT of the OAF protocol 
is identical to the DMT of the Non-Orthogonal Amplify and 
Forward (NAF) protocol. 

Two variants of the NSDF protocol are considered: flxed-NSDF 
and variable-NSDF protocol. In the variable-NSDF protocol, the 
fraction of time duration for which the source alone transmits is 
allowed to vary with the rate of communication. Among the class 
of static amplify-and-forward and decode-and-forward protocols, 
the variable-NSDF protocol is shown to have the best known 
DMT for any number of relays apart from the two-relay case. 
When there are two relays, the variable-NSDF protocol is shown 
to improve on the DMT of the best previously-known protocol for 
higher values of the multiplexing gain. Our results also establish 
that the flxed-NSDF protocol has a better DMT than the NAF 
protocol for any number of relays. 

Finally, we present a DMT optimal code construction for the 
NAF protocol. 

Index Terms — cooperative diversity, distributed space-time 
code, orthogonal amplify and forward, non orthogonal am- 
plify and forward, selection decode and forward, diversity- 
multiplexing gain tradeoff, space-time codes, cyclic division 
algebra codes. 

I. Introduction 

Cooperative relay communication is a promising means of 
wireless communication in which cooperation is used to create 
a virtual transmit array between the source and the destination, 
thereby providing the much-needed diversity to combat the 
fading channel. 

Consider a communication system as shown in Fig. Q] in 
which there are n + 1 nodes that cooperate in the communi- 
cation between source node S and destination node D. The 
remaining (n — 1) nodes thus act as relays. 
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Fig. 1. Cooperative Relaying in networks. 



A. Assumptions 

We follow the literature in making the assumptions listed 
below concerning the channel. Our description is in terms of 
the equivalent complex-baseband, discrete-time channel. 

• All nodes have a single transmit and single receive 
antenna and are assumed to transmit synchronously. 

• The T channel uses over which communication takes 
place is short enough to invoke the quasi-static assump- 
tion, i.e., the channel fading coefficients {g i: hf\ are fixed 
for the duration of the communication, but vary randomly 
from one block of T channel uses to the next. Here, g\ 
denotes the fading coefficient for the channel between S 
and D, gt for the channel between S and the ith relay 
node Ri, and hi for the channel between R4 and D, where 
2 < i < n. 

• All channels are assumed to be Rayleigh fading and all 
fade coefficients are i.i.d., circularly symmetric complex 
gaussian CA/"(0, 1) random variables. 

• We assume half-duplex operation at each node, i.e., at 
any given instant a node can either transmit or receive, 
but not do both. 

• The noise vector at the receivers is assumed to be com- 
prised of i.i.d., circularly symmetric complex gaussian 
CA/"(0, a 2 ) random variables. 

• The destination must know all the fading coefficients in 
the case of the AF protocols and it needs to know only the 
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source-destination and relay-destination channel fading 
coefficients in the case of the DF protocols. We assume 
that a relay knows only the corresponding source-relay 
channel fading coefficient, i.e., relay Rj will know only 
gj in the case of the DF protocols. 

B. Channel Model 

We propose several cooperative communication protocols 
in this paper. All the protocols considered here involve two- 
phase communication. In the first phase, lasting for p channel 
uses, the source S broadcasts to the relays and the destination. 
In the second phase, lasting for q channel uses, the relays 
communicate with the destination. A protocol is said to be 
non-orthogonal or orthogonal depending on whether the source 
continues to transmit (to the destination) in the second phase 
or otherwise. It is said to be a decode-and-forward (DF) or 
amplify-and-forward (AF) protocol depending upon whether 
the relays are required to decode the received message or not. 

As we shall see, one can associate with every protocol, a 
channel model of the form 

y = Hx + w, ( l ) 

where y corresponds to the received signal, w is the noise 
vector, H is the virtual channel matrix and x is the vector 
transmitted by the source in the case of the AF protocols, 
and is the compound vector formed by concatenating the 
transmissions of the source and the participating relays in case 
of the DF protocols. 

C. Diversity-Multiplexing Gain Tradeoff 

The average signal-to-noise ratio of a channel will be 
denoted by p. In the channel model in (03, the communication 
has happened over m = p + q channel uses. 

The probability of outage for the channel in (Q]i is defined 

as 

P aat {R)= mf Pv(I(x;y\H) < mR), 

S x > 0, Tr(£ x ) < p£ - 

where R is the average rate of communication. We have 
imposed an average energy constraint on the transmitted signal 
x by upper bounding the trace of its co variance matrix S x . As 
will be shown later, by normalizing the noise vector, we can 
regard £ as p. 

In this paper, we use the symbol = to denote exponential 
equality, i.e., the expression 

l im ^M = 6 

p^oo log/9 

is denoted by f(p) = p b and >, < are similarly defined. 

We define d out (r) to be the negative p exponent of 
P out (r log p), where R = rlogp, i.e., 

Pout(r) =p-**M. 

A MIMO coding scheme {C(p)} for the channel model in 
(Q~|i is said to achieve a spatial multiplexing gain r if 

r R(p) 

lim = r 

p-+oo logp 



where R(p) is the rate of the code C(p). The coding scheme 
{C(p)} is said to achieve a diversity gain d(r) if 

lim X ° MP) = -d(r) , 

p—>oo log p 

where P e (p) is the average error probability of the code C(p) 
under maximum likelihood decoding. Thus 

P e {p) = p- d ^. 

It can be shown from an application of Fano's inequality (see 
[2]) that 

d(r) < d oat (r). 

The diversity-multiplexing gain tradeoff (DMT) introduced 
by Zheng and Tse in [2], provides a means of evaluating and 
comparing the various proposed protocols. In this paper, we 
will by abuse of terminology refer to the outage exponent 
d ut(r) as the DMT of the corresponding channel even though 
strictly speaking it is the DMT only if one can identify a 
code whose diversity gain d(r) is equal to d out (r). No loss 
of accuracy is entailed in this because, as we shall see, for 
every protocol discussed in this paper we are able to identify 
a correspondingly optimal coding scheme. For this reason, we 
often will write d(r) in place of the more accurate, but more 
cumbersome notation d out (r). 

For each of the protocols described in the paper, the relation 
between received and transmitted vectors can alternately be 
written in the form 

Y = HX + W , 

where H is a row matrix whose components represent the 
fading coefficients a MISO channel and where X is a code 
matrix, constrained by the particular protocol, to have a 
prescribed format. From this, it follows that the probability of 
outage of the channel model in ([T]i is at least as large as that 
of the corresponding MISO channel. Since the MISO channel 
is known to have outage exponent given by 

d(r) = n(l — r) , 

when (n — 1) is the total number of relays, the quantity d(r) 
is referred to in the literature as the transmit diversity bound 
[4]. 

D. Prior Work in Cooperative Communication Protocols 

The concept of user cooperative diversity was introduced in 
[12], [13]. Cooperative diversity protocols were first discussed 
in [4], where the authors develop and analyze the Orthogonal 
Amplify and Forward (OAF) protocol and the Selection De- 
code and Forward (SDF) protocol for the case of a single 
relay. In [6], the SDF protocol is analyzed for an arbitrary 
number of relays, where the authors give upper and lower 
bounds on the DMT of the protocol. In these protocols, the 
relays and the source node participate for equal time instants 
and the maximum multiplexing gain r that could be achieved 
was 0.5. 

In [5], Azarian et al. analyze the class of Non Orthogonal 
amplify and Forward (NAF) protocols, introduced earlier by 
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Nabar et al. in [11]. In [5], the authors establish the improved 
DMT of the NAF protocol in comparison to the class of OAF 
protocols considered in [4]. The authors also introduce the 
Dynamic Decode and Forward (DDF) protocol wherein the 
time for which the relays listen to the source depends on the 
source-relay channel gain. They show that the DMT of the 
DDF protocol achieves the transmit diversity bound for r < 
0.5, beyond which the DMT falls below the transmit diversity 
bound. 

Jing and Hassibi [7] consider cooperative communication 
protocols where the relay nodes apply a linear transformation 
to the received signal. The authors consider the case when 
both the source and the relays transmit for an equal number 
of channel uses and the linear transformation applied by the 
relays are restricted to the class of unitary matrices. 

Yang and Belfiore consider a class of protocols called 
Slotted Amplify And Forward (SAF) protocols in [10], and 
show that these improve upon the performance of the NAF 
protocol [5] for the case of two relays. The authors also 
provide an upper bound on the DMT of the SAF protocol with 
any number of slots, and show that this upper bound tends 
towards the transmit diversity bound as the number of slots 
increases. Under the assumption of relay isolation and relay 
ordering, the naive SAF scheme proposed in [10] is shown to 
achieve the SAF protocol upper bound. 

Yuksel and Erkip in [18] have considered the DMT of 
the DF and compress-and-forward (CF) protocols. They show 
that the CF protocol achieves the transmit diversity bound 
for the case of a single relay. We note however, that in the 
CF protocol, the relays are assumed to know all the fading 
coefficients in the system. 

Recent work directed towards the construction of explicit 
coding schemes which are DMT optimal for cooperative 
protocols can be found in [9] and [10]. 

In the present paper, neither relay ordering nor relay isola- 
tion is assumed. 

E. Results 

The results presented in this paper pertain to four protocols. 

1) Orthogonal Amplify and Forward: This protocol was 
introduced by Laneman et al. [4], who analyze this protocol for 
the case of a single relay. As the name of the protocol suggests, 
in this two-phase protocol, the source broadcasts for p time 
slots, followed by a relaying phase, lasting for q time slots, 
where the relays amplify and forward the signals received by 
them. 

Our version of the protocol is slightly more general than that 
in [4] since we permit the relays to operate on the received 
signal using a linear transformatiorQ and we allow the source 
and the relays to transmit for unequal time slots. For this more 
general version of the protocol, we are able to determine the 
best possible DMT as well as construct, in a simple way, DMT 
optimal codes that incur minimum delay. Our results establish 
that the DMT of the OAF protocol matches the best possible 
DMT of the NAF protocol [5]. 

In this sense, a more appropriate term for this protocol we consider here 
might be "linearly transform and forward". However we shall use the more 
generic term OAF here. 



2) Selection Decode and Forward: This class of protocols 
was introduced by Laneman and Wornell [6]. In this protocol, 
in the first phase occupying p channel uses, the source broad- 
casts to the destination and the relays. In the second phase, all 
the relays which are not in outage decode the source message, 
separately encode and transmit to the destination over the 
course of q channel uses. Once again, we consider the case 
when p q. 

We consider two versions of this protocol in the present 
paper: 

• Non-Orthogonal Selection Decode and Forward (NSDF) 
In this protocol, the source continues to transmit during 
the second phase. At this point, it is convenient to 
distinguish between two variants of this protocol. 

- Variable-NSDF In this variant, in order to obtain the 
best DMT possible, we allow p and q to vary with 
the multiplexing gain r. Note that since, p and q 
are not a function of the channel fading coefficients, 
this protocol still falls within the category of static 
protocols. 

* Among the class of static DF and AF protocols, 
the variable-NSDF protocol has the best-known 
DMT for any number of relays, except for the 
case of two relays. 

* The DMT of the variable-NSDF protocol for the 
case of two relays is better than the tradeoff of 
the SAF protocol [10] for r > 0.6 (refer Fig. |3). 

- Fixed-NSDF Here p and q are fixed and independent 
of r. The DMT for this variant of the NSDF protocol 
is also presented here for every pair p,q, p> q. 

* For values of the ratio k = £ in the range 1 < 
n < ^zj, the fixed-NSDF protocol dominates the 
NAF protocol for < r < beyond which 
they both have the same DMT. 

• Orthogonal Selection Decode and Forward ( OSDF) 

In this protocol, the source remains silent in the second 
phase. As in the case of the NSDF protocol, there are 
two variants of this protocol, fixed and variable. Under 
the variable-OSDF protocol, we allow p and q to vary 
with r in order to compute the best DMT. We state the 
DMT for both the fixed and variable-OSDF protocols. 
The DMT of the variable-OSDF protocol, for the case of 
two relays, improves on the tradeoff of the SAF protocol 
[10] for r > §. The DMT of the variable-NSDF protocol 
is, however, better than the DMT of the variable-OSDF 
protocol for all r for any number of relays. 
All the OAF, NSDF and OSDF protocols considered in this 
paper are static protocols. In Fig. [2] and Fig. [3] we have shown 
the optimal DMT of the OAF, NSDF and OSDF protocols for 
the case of one and two relays respectively. In the figures, we 
have also shown the DMT of the NAF and SAF protocols for 
the sake of comparison. Note that the NAF and SAF protocols 
have the same DMT for the case of a single relay. 

As our final result, we present a code that achieves the DMT 
of the NAF protocol considered by Azarian et. al. [5]. This 
result, first presented in [17], is included here for the sake of 
completeness. 





Fig. 3. Optimal DMT for two relay cooperative communication protocols. 
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F. Organization of the Paper 

In Section [LI] we discuss the class of OAF protocols and 
compute their DMT. Section [HI] contains the description and 
DMT analysis of the NSDF protocols. We state the DMT of the 
OSDF protocols in Section [IV] We also present constructions 
for DMT optimal codes, based on cyclic division algebras, for 
the OAF, NSDF and OSDF protocols in the corresponding 
sections. Finally, we provide a code construction which is 
DMT optimal for the NAF Protocol [5] in Section [V] 

Notation: The norm of a vector and the Frobenius norm 
of a matrix are denoted by j | . 1 1 and \\.\\f respectively. | . | will 
denote the determinant of a matrix as well as the modulus of 
a complex scalar. 

II. Orthogonal Amplify and Forward Protocol 

In this protocol, the source £ transmits a signal to the relays 
{Rj} and to the destination D for p channel uses. Over the 
next q channel uses, the relays transmit a linear transformation 
of the received signal, while the source remains silent. 

As it turns out, for -7- = „ - ' , we get the best DMT for 

' P+q 2n— 1' & 

all r in the range < r < 1/2. For r > 1/2, we obtain the 
best DMT when there is no cooperation from the relays and 
the source continuously transmits to the destination. 

Hence, we get the DMT of the OAF protocol by choosing 

1 when r > 1/2. 



p 

p+q 



2n-l 



when r < 1/2 and 

— ' p+q 



A. OAF Channel Model 

Based on the above signalling protocol, we have the fol- 
lowing model for the received signal. 



r : 



H2 



gix + w ± 

gjx + Vj , j = 2,3, ... ,71 



J=2 



J=2 



x + ^hjAjVj 

j=2 



W 9 



where x is the signal transmitted by the source, [y*,?/^]' * s 
the signal received by the destination and r_j is the signal 
received by the jth relay. {Aj} are (q x p) complex matrices 
that represent the linear transformation at the relay nodes. The 
vectors {Vj}™—2 and {w 1; w 2 }' represent additive noise seen 
at the relay nodes and the destination respectively. The above 
channels can be rewritten in matrix format as 



V = Hx + n, 



(2) 



The covariance matrices of the noise and signal vector are 
denoted by 



5j« 



i 

v- 





^■2 T 



°llq 



En 
J=2 



■AjA] 
1 j j 



and 



T, x := E(xx}), 



where a% . , tr^ denote the variances of the corresponding noise 
vectors. We impose the energy constraint, 

Tr^) < P S , 

where £ denotes the average energy available for transmission 
of a source symbol. The SNR p is defined as 



£ 

P '■= — 



(3) 



We assume the ratio of the noise variances {a^. } and to be 
a constant independent of p. The average energy of the signal 
transmitted by the j relay is given by, 



E{\\g J A J x + A.vA] 2 } 



= \gtf EiWAjxf} +Tr(A j E{v j J j }A]) 
< | ft |VTr{^4} + ^Tr{^4} 



< ^{p£\g 3 \' + a z w ) 

where is the squared Frobenius norm of the relay matrices 
{Aj}, i.e., 



(4) 



\\A^ F = Tt{AiA)) := a). 
We have used the inequality 

£2 < p£l P 

in coming up with the above bound. 

We impose the constraint that the average energy transmit- 
ted by a relay satisfy 

EiWgjAjX + A^W 2 } <£, 

and this is achieved by requiring that 



p£\gj 



(5) 



As in [5], the SNR exponent of a 2 can be made zero. In 
practice, we can choose aj to be a suitable constant so that 
the probability of the event described by © is high. 



where 



H = 



E"=2 hjAjVj 

9i!p 
X^=2 9jhjAj 



Hi 



B. Upper bound on the DMT of the class of OAF protocols 

We first state a useful lemma concerning non-negative 
definite matrices. 

Lemma 1: Let {gj, hj, Aj} be as defined above. Let 



B =^2,9jhjAj. 

3=2 
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Then, 



BB* < (n- 1)^7^5 

3=2 



as non-negative definite matrices, where jj = \gjhj\ 2 . 

Proof: Please see Appendix U ■ 

Next, we establish the upper bound on the DMT of the class 
of OAF protocols. 

Theorem 2: (General OAF Upper Bound) Consider the col- 
lection of OAF protocols described above (different protocols 
can be obtained by varying p, q and {A,} for a given n). 
Then, regardless of the choice of the transformation matrices 
{Aj}, the DMT of any protocol satisfies the upper bounds 
given below. 

If „ > o„" -i , where m = p + q, 



d(r) 



2n-l ' 



< 



n 1 



(n-l)mr 
nq 



p / ^ mr 

p-q \ p 

(1-0 , 



< r < 

— — m 
<r < i 

m — — 2 

1 < r < 1 



(6) 



If x- < 



2n-X 



then 



d(r) 



< 



n(l 
(1 



P 



) , < r < 



("-!) 



") . 



< r < 1 



(7) 



In deriving these bounds for the protocols, cooperative 
relaying is avoided whenever it is advantageous to do so. 

Also, the highest value of the upper bound on the DMT 
occurs for the choice X. — _» t j n this case, we get 



d(r) 



< 



2n 

1 _(2n^)r 
n 

(l-r) , 



0<r<i 
i < r < 1 



(8) 



Proof: As in [2], the input distribution can be assumed 
to be circularly symmetric gaussian without loss of generality. 
The maximum mutual information between x and y, condi- 
tioned on the knowledge of H at the receiver, is given by 



max I(x: y I H) 

E x > 0, Tr(E x ) < p£ - 

max log I /„ 

S x > 0, Tr(£ x ) < p£ 



Arguing as in [2], in the scale of interest we get 

2max = log llm+SHH^- 1 \ 

It turns out that exponential equality as described above is 
sufficient to determine the outage probability in the scale of 
interest. Let 



J 



P9i 



pg\B Iq + pBB^C- 1 



where C = I g + EL 2 W^AjA], 71 = |. 9 i| 2 and B is as 
defined in Lemma Q] 



Then upon row reduction of J, we obtain 

I J I - 1(1 + Pli)I P I • I Iq+ —^BB^C" 1 I 

1 + P7i 

= (I + ptOHC- 1 I \C+—^- Bfltj (9 ) 

1 + P7l 



By applying Lemma [T] and using ®, we get 



1^1 < (l + P7i) P - (i + E 

3=2 



1^1 



P 



1 + P7i 



7j 



< (l + P7i) P_9 ( 1 + P7i+P^7j 

J=2 



As in [5], we define 

7i 

For j = 2, 3, . . . , n, 



j~ u , u>0 
0, u < 



7j 



(10) 



(ID 



and 



v = min{«j} j=2 . 



This gives us 

J J J < /9 (p-?)(l-i') + +gmax{(l-ti),(l-t))} + ^ 

The probability of outage, for the channel in (O, is defined as 

Pout(R)= inf Pr(I{x;y\H) < mR). 

S x > 0, Tr(S x ) < p£ ~ 

In the scale of interest, the above expression reduces to (see 
[2]), 

Pout(R = r\ogp) = Pr( log I J I < rm log p ) 
and is lower bounded by 

Pr ((p — q)(l — u) + + q max{(l — u), (1 — v)} + < rm) . 
Let the negative exponent of P out (R) be d(r), 

log Pom (r log p) 



d(r) 



lim 

p — >oo 



logp 



Then, 

d(r) < 



inf u -\- {tl — l)v. 

(p— g)(l — + max{(l — — < mr 



It is clear that it is enough to consider u, v < 1. Hereafter, 

we will consider u and w to lie in the range < u, v < 1. 
Therefore, 

< inf u + (n — l)v. 

(p—q)u+qmin{u,v} > p—mr 

By solving the above optimization problem, we get the 
statement of the theorem. Please see Appendix [TT] for the 
details. ■ 
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C. Specific protocol achieving the DMT upper bound for p = 

n, q = (n - 1) 

Within the class of OAF protocols, there are different pro- 
tocols corresponding to various choices of p, q and {A,}™ =2 . 
As seen in Section ITl-BI for a given number (n— 1) of relays, 
the upper bound on d(r) is maximized when = 2 n-i • 

We have used the inequality 



9jhjAj 




uJ'j I 2 AiA\ (12) 



to derive an upper bound on the DMT of the class of OAF 
protocols. Equality will occur in dT2l > if 



AjA\ = , for all j ^ k, 



(13) 



i.e., if the row spaces of the matrices {Aj}™ =2 ^ pairwise 
orthogonal. 

The arguments presented above serve as a motivation for the 
particular choice of the matrices {Aj} outlined next. With this 
specific choice of {A,}, satisfying the constraints in (1131 1. we 
can achieve the upper bound on the DMT given in Theorem 

m 

Theorem 3: Consider a specific OAF protocol, as described 
above, with parameters p = n and q = n — 1. Choose the 
(n-l)xn matrices {Aj} as follows: 



Mk,i) 



o 



k = j -1,1 = j 
elsewhere 



(14) 



i.e., the (j — entry of Aj is equal to aj and remaining 

entries are 0. The DMT of this protocol is equal to the highest 
upper bound of the class of OAF protocols, 



d(r) 



(2n- l)r 
1- r 



< r < i 
h < r < 1 



Proof: With the above choice of {Aj} we get, 
BB^ 



By substituting for BB^ in equation (O, and setting p = n 
and q = n — 1, we get, 



J 



(1 + phT- 

1-L_£ 



1+P71 



72 



+ 1+P71 7 " 



(1 + P7i) + Pli + PI]) 

j=2 



J=2 



-min{ u,Vj })"*" 



As in Theorem |2] the outage probability is given by 

■Pout (flog p) 

= Pr(log \J\< (2n-l)r\ogp) 
( 

= Pr (1 - u) + + (1 - min{u, v 3 }) + < (2n - l)r 

Let the negative exponent of P out (rlog/o) in the scale of 
interest be d(r). Then, we can consider u and to lie in 
the range < u, Vj < 1. We get 



d(r) = inf u 

u+^7J" =2 min{ii,i)-/}>n— (2n— l)r 

By solving the optimization problem, we get 

d(r) = n - (2n - l)r , < r < 



J=2 



2n- 1 



(15) 



For r > i, we allow the source to transmit to the destination 
continuously, thereby achieving the tradeoff mentioned in the 
statement of the theorem. ■ 

Since 2n— 1) = 1, the smallest value of delay parameter 
m = (p + q) satisfying the condition = 2n "_ 1 corresponds 
to the choice p = n, q = (n — 1). Hence, the above protocol 
has minimum possible delay required to achieve the best DMT. 

Example 1: Let the number of relays be 2. Therefore, n = 
3. We choose q = 2, p = 3, and 



OL 2 





A, 





a 3 



For these parameters, the DMT of the OAF protocol is 

d(r) = 
and is shown in Fig[3] 



3 - 5r 
1 -r 



. 0<r<i 
. 5<r<l 



D. DMT Optimal Codes for the OAF Protocol 

Our code construction is based on cyclic division algebras 
(CDA) and we begin with a brief introduction to these alge- 
braic objects. 

1) Division Algebras: Division algebras are rings with 
identity in which every nonzero element has a multiplicative 
inverse. The center F of any division algebra D, i.e., the subset 
comprising of all elements in D that commute with every 
element of D, is a field. The division algebra is a vector 
space over the center F of dimension n 2 for some integer 
n. A field L such that FcLcfl and such that no subfield 
of D contains L is called a maximal subfield of D (Fig. 0). 
Every division algebra is also a vector space over a maximal 
subfield and the dimension of this vector space is the same for 
all maximal subfields and equal to n. This common dimension 
n is known as the index of the division algebra. 
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D 



Division Algebra 



Maximal Subfield 



F Centre 
Fig. 4. Structure of a division algebra. 



2) Cyclic Division Algebras: Our interest is in CDA, i.e., 
division algebras in which the center F and a maximum 
subfield L are such that L/F is a finite cyclic Galois exten- 
sion. CDAs have a simple characterization that aids in their 
construction, see [15, Proposition 11], or [14, Theorem 1]. 

Let F, L be number fields, with L a finite, cyclic Galois 
extension of F of degree n. Let a denote the generator of the 
Galois group Gal(L/F). Let z be an indeterminate satisfying 



Iz 



za(£) V feL and 



for some non-norm element 7 6 F*, by which we mean some 
element 7 having the property that the smallest positive integer 
t for which 7* is the relative norm N^^ ¥ (u) of some element 
u in L*, is n. Then, a CDA D(L/¥, a, 7) with index n, center 
F and maximal subfield L is the set of all elements of the form 



\d under this basis has the matrix representation, 



£ 7^(4-i) 7^ 2 (C- 2 ) 



7a 











(17) 

known as the left regular representation of d. 

A set of such matrices, obtained by choosing a finite subset 
of elements in D constitutes the CDA-based ST code X. In 
[1], the authors have constructed CDA-based ST code for all 
values of n. For all the codes constructed in [1], the underlying 
constellation is QAM and the center of the division algebra is 
F = Q[i}. 

4) DMT optimal CDA-based ST codes for the OAF protocol 
with parameters p = n,q = (n — 1): In this subsection, we 
provide an explicit construction of a code, based on CDAs, 
for the OAF protocol for any number of relays, and prove 
the DMT optimality of the code. If the number of relays is 
n — 1, we choose p = n and q = (n — 1) since this choice of 
parameters has the best DMT (see Theorem O. 

The channel model for the OAF protocol in (ffji can be 
rewritten as 



[ id & ) 



[ 91 92112 











(18) 



n-l 

E 

i=0 



z l L 



L G L. 



(16) 



Moreover, it is known that every CDA has this structure. It can 
be verified that D is a right vector space (i.e., scalars multiply 
vectors from the right) over the maximal subfield L. 

3) Space-Time Codes from Cyclic Division Algebras: A 
space-time (ST) code X can be associated to D by selecting 
the set of matrices corresponding to the matrix representation 
of elements of a finite subset of D. Note that since these ma- 
trices are all square matrices, the resultant ST code necessarily 
has T = n t . 

The matrix corresponding to an element d E D corresponds 
to left multiplication by the element d in the division algebra. 
Let Ad denote this operation, Ad : D — > D, defined by 

A d (e) = de, V e G D. 



where 



For M even, let „4 QAM denote the M 2 -QAM constellation 
given by 



A QAM = {a + ib I \a\,\b\ < M - 1, a,b odd} . 



(19) 



Consider a CDA having center F = Q(i) and maximum 
subfield L that is a degree-?i cyclic Galois extension L/F of F. 
Let a be the generator of the cyclic Galois group Gal(L/F). 
Let Or and Ol denote the ring of algebraic integers in F and 
L respectively. It is known that Of = Z[i]. Let {/3i, . . . , /?„} 
be an integral basis for Ol/Of- Let D(L/F, a, 7) denote the 
associated CDA. 
Let 



G -4.qam(/3i, ' ' ' 5 (3n) 



(20) 



It can be verified that Ad is a L-linear transformation of D. 
From dl6l . a natural choice of basis for the right-vector space 
D over L is {1, z, z 2 , . . . , z" -1 }. A typical element in the 
division algebra D is d = £q + z£i + ■ ■ ■ + z" _1 £„_ 1 , where 
the £i G L. By considering the effect of multiplying dx 1, dx z, 
. . . , dx z n_1 . one can show that the L-linear transformation 



where 



^4qam(/5i , ■ ■ ■ j 0n) 



G *<4q A 



Then, let the signal transmitted by the source in the first 



9 



p = n channel uses be given b>0 



x = [to cr{£o) ... <r n -\l a ) Y (21) 
where £ G Ajam(A: • • • , /?«)• 

We select {A,} as specified in ( fl"4l i. Without loss of 
generality (insofar as DMT is concerned), for simplicity we 
set 

a.j = 1 , 2 < j < n. 

As we saw earlier in Theorem[3] this choice of {Aj} achieves 
the upper bound on the DMT for the OAF protocol. Hence, 
the code matrices will be as shown in d22l . 

The performance of this code is no worse than that obtained 
by deleting the columns 2 to n — 1 of all the code matrices X. 
Hence, from here on, we will work with the column-deleted 
code matrix: 



X 



L (4)) 



where £ G A>am(/3i, • • • , /?«)■ 



Theorem 4: The DMT of the above code is 



(23) 



d(r) 



(2n - l)r , < r < 



2n- 1 



Proof: The mentioned tradeoff is the optimal DMT for 
the OAF protocol in the range < r < |. 

For choice of the Aj matrices given by ( TBI ), the covariance 
matrix of the noise vector n is shown in ( f24b . By normalizing 
the noise variances to unity and noting that 1 + \ hj\ 2 = 1, the 
covariance matrix becomes identity in the scale of interest. 
Hence, we can consider the noise to be white in the scale 
of interest. Also, we impose an energy constraint on the 
codewords and replace X with OX, where 9 is chosen to 
ensure that 

\\6Xf F < (2»-l)p. 

After deleting columns 2 to n — 1 of the ST code for the 
OAF protocol, the resultant channel model can be rewritten as 



y = 



91 



g n h n 

where £ G A QAM (Pi, 



to 



l (£o) 



+ n 



It can be shown that, in the scale of interest, the above channel 
is equivalent to a parallel channel. Hence, we show that the 
chosen code is an optimal code for the parallel channel. 



To support a data rate of R p 



, log(p) on the parallel 



channel, we need, M = p » . The energy requirement forces, 



2 £q need not be drawn from a maximal subfield L of a division algebra. 
It is enough if L is an algebraic extension of Q[i] of degree n and £q is as 
defined in (21) , 



Now, the product of the squared norms of the normalized 
difference code matrices (obtained by scaling the code vectors 
with -^=) is given by 



e \*-\Oa(£o)f---\0<T n - 1 (to)\ 2 



rfln n ~ 1 

^ni^°)i 



2 > P '_ p 



i=0 



Therefore, from [3, Theorem 5.1] the code is DMT optimal 
for the parallel channel. It can be shown that the DMT of the 
parallel channel is (see [3]) 



dp {j"p ) ^ ^p •} 



< r p < n. 



(25) 



Let the rate of the original code for the OAF protocol be 
R = r log p. Then, the size of the DMT optimal code book 
for the parallel channel is p Tp . This corresponds to a code 
book for the OAF channel of size ( o( 2 ™~ 1 ) r . Since the two 
code books are the same, it follows that r p = (2n — l)r. By 
substituting for r p in (l25l l, we obtain a lower bound on the 
DMT of the code for the OAF protocol. The lower bound 
occurs because dropping some columns of the code matrices 
could conceivably decrease the probability of error. Therefore, 

d(r) > d p ((2n-l)r) 

n 

> n - (2ra- l)r , < r < 



2n-l 

The above bound equals the DMT of the OAF protocol in 
(15[ . Therefore, the DMT of the code is as mentioned in the 
theorem. 

For r > i, the source will transmit continuously to the 
destination and the relays will not participate. It follows that 
the code is DMT optimal for the class of OAF protocols 
considered here. ■ 

Remarks: We make the following remarks on the class of 
OAF protocols and the proposed DMT optimal code: 

1) Among the class of OAF protocols, the best DMT is 



achieved when -2- 



2n-l 



Since n and 2n — 1 are 



relatively prime, the proposed code has the minimum 
possible delay of 2n — 1 with the parameters p = n and 
q = n — 1. 

2) For p = n and q = n—1, the DMT of the OAF protocol, 
with the choice of {Aj} mentioned in Theorem [3] 
coincides with the DMT of the NAF protocol [5]. 

3) When each node in the system has only one transmit 
and one receive antenna, the DMT optimal code for the 
NAF protocol proposed in [9] has delay 4(n— 1) which 
is larger than the delay for the DMT optimal code for 
the OAF protocol. 

III. Non-Orthogonal Selection Decode and 
Forward Protocol 

In this section, we consider a non-orthogonal selection 
decode and forward (NSDF) protocol in which the source 
transmits a signal to the destination and the relays for p 
channel uses in the first phase. All the relays, which are not 
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4 <t(40 




o 



X = 







(4) 



where 4 e Ajam(/3i, ■ • • ,/?«)• 



(22) 



cC + (J. 



?L? 1 Vn 



(24) 



in outageQ, will decode the source message. In the second 
phase, the relays will separately encode and transmit a vector 
of length q. The source continues to transmit to the destination 
in the second phase. We only consider the case when p > q. 
To compute the best possible DMT, we allow p and q to vary 
with the multiplexing gain r and choose the value of k = | 
which maximizes the DMT for a given r. This version of the 
protocol will be called the variable-NSDF protocol. We have 
also computed the DMT of the fixed-NSDF protocol, wherein 
the ratio k = | is fixed for all r. We have constructed CDA 
based codes which achieve the DMT of the variable and fixed 
NSDF protocols. 



A. DMT of NSDF Protocol 

Theorem 5: The DMT of the variable-NSDF protocol is 
given by 



d(r) 



/-, (n-l)(«„+l) 

n 1 — - — i 

n 

(n — r)(l — r) 
(n-2)r+l ' 



< r < — K-r 

— — re n +l 

< r < 1 



K„ + l 



(26) 



where n n = v 2(n _ 1} . 

In deriving the above DMT, we have allowed p and q to 
vary with the multiplexing gain r. The source and the relays 
choose a code corresponding to each p and q. For the case 
p > q, we select the value of (p, q) which maximizes the 
DMT for a given r. Suppose p = nq, then the optimal value 
of k is given by 

f Kn, < r < -At 



l+(n-2)r 1 
(n-l)(l-r) ' k„ + 1 



< r < 1 



For a fixed choice k 
protocol is given by: 

if 1 < k < K n , 



the DMT of the fixed-NSDF 



d(r) = (n-1) ( 1 - — ) +(l-r) , < r < 1 , (27) 



3 We say that a relay is not in outage if the corresponding source-relay 
channel is not in outage. 



else if k > k,- 



d(r) 



1 - 



m(n — 1) ( 



Tig 

f (1-r) , 

_ (ra— l)m+p 



< i 
— m 

np—m 



i 



- r 



< r 
< r < 

— — (n— 2)m+p 
np-m < r < ii 



(n— 2)m+p 

i < r < 1 



(28) 



Proof: Consider the system model described in Section U 
Let xi and x 1 be the signals transmitted by the source in the 
first and second phase respectively. Let (fc — 1) relays, where 
1 < k < n, participate in the second phase and let {xj}j =2 
be the signal transmitted by them. Only the relays that are not 
in outage in the first phase shall participate in the cooperative 
protocol in the second phase. 

Let Ek denote the event when any (fc — 1) relays 
participate in the cooperative (second) phase. The events 
{Ei , E 2l ...,£"„} are disjoint and their probabilities sum up to 
1. First, we will calculate the outage probability of the channel 
when the event Ek occurs. 

1) Outage Probability conditioned on Ek-' Consider the 
case when the event Ek, 2 < fc < n, occurs. The case fc = 1 
will be dealt with separately. Let the signals received by the 
destination in the two phases be yi and y 2 , where 



gixi + wi 



and 



2/2 



gixi 



k 

i=2 



W2 



with wi, W2 denoting the noise added at the destination in the 
respective phases. We impose an energy constraint by choosing 
£ to be the average energy available for transmission of a 
symbol at either the source or a relay. Let cr^ be the variance 
of the noise added at the destination and SNR be as defined 
in ([3]). We shall normalize the noise variances to unity. Hence, 
we can regard £ as the SNR p. 

The channel model for the NSDF protocol is given in (f29b - 
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Ml 




9il, 


. 2/2 _ 










9llq 



£i 
■11 

X2 



Wl 
W2 



(29) 



Let, 



and 



U. 

H k = 



Ml 
Ml _ 
t '* 

X\ x 1 



Wl 

w 2 



giip 







x 2 











h 2 I q ■ 

Then, ( 1291 can be rewritten as 

y = H k x + w. 

The input vector can be assumed to be circularly symmetric 
complex gaussian. The maximum mutual information trans- 
ferred between the source and the destination, conditioned on 
the knowledge of H k at the destination, is given by 

I{mx\H k )=\og\I m +H k Z x Hl | . 

Arguing as in case of the OAF protocol, in the scale of interest, 
the maximum mutual information between x and y is given 
by 



log | I m + P H k Hl | 

k 

\og(l + p\g 1 \ 2 ) P (l+p\g 1 \ 2 + P Y / \h ] \ 2 ) q 

i=2 

log// (1 ~" )++9max{(1 ~' ii) ' (1 ~' u)}+ 



where, 



l.9i I = 



N 2 



u > 
u < 

Vj >0 
Vj <0 



and 



l i v 3}j=2- 



The outage probability of the channel in 
of interest, is given by, 

-Pout (r log p) = Pr(X max < mr log p) 

so that, 



, in the scale 



:= p~ dk( - r ^ 



d k (r) = inf u+(k-l)v. (30) 

p(l — u)+ -\-q max{(l-n),(l-u)}+ <mr 

It is clear that it is enough to consider < u, v < 1. Hence, 
the above infimum must be calculated over the region 

pu + gmin{u, v} > m(l — r) , < u, v < 1. 



It is clear that for r > 1, dk{r) = 0. Therefore, it is enough 
to consider r < 1. We consider two separate cases to evaluate 
dk{r). We set p > q in the optimization that follows. 
Case I: min{u, v} = u 
We get, 

u > (1 — r) and v > u. 



Substituting in d30t we get, 

d k (r) = fc(l-r) , < r < 1. 

Case II: min{u, u} = v 
We have, 

pu + qv > m{l — r). 

As in the case of the OAF protocol, we solve the optimization 
problem to get, 

f t_ (fc-l)m r n < r < -2- 



Now, we handle the case when event E\ occurs, i.e., k = 1. 
The channel model in this case is given by 



Ml 

</2 



9ii P 
o 





9l!q 







+ 










U> 2 



(32) 



= H\x + w. 

The probability of outage of the above channel is given by 

Pout(rlogp) := p- d ^ =p-^ + . (33) 

2) Probability of the set of Participating Relays: In this 
subsection, we will compute the probability of the event E k . 
The probability of E k is the product of probabilities of two 
events: 

Pi(E k ) = Pr((n — k) relays are in outage) • 

Pr((fc — 1) relays are not in outage). 

We shall evaluate the probabilities of the two events mentioned 
above separately. Also, we say that a relay is participating in 
the second phase if the corresponding source-relay channel is 
not in outage. 

The signal received by the jth relay Rj, where 2 < j < n, 
in the first phase is given by 

Tj_ = 9jXi + Vj_ , 

where Vj is the noise vector. The maximum mutual informa- 
tion between x\ and Vj, conditioned on the knowledge of gj 
at the relay, is 

J max = max I(rj-,xi\gj) 

S X1 > 0, Tr(S xl ) < p£ — 

V 



m 



log(l + p|&f 
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Now, the probability that Rj is in outage is 
Pr(i? 3 is in outage) = Pr ^— log(l + p\gj\ 2 ) < rlogp 



where 



\9j\ 



rr^i- Uj y< — 



P~ U] , Uj > 
, Uj < 



Since the fading coefficients corresponding to different source- 
relay channels are independent of each other, we have 

Pr((n - k) relays are in outage) = p - ^ - *^ 1-2 ?) . 

The probability of a relay participating in the second phase 
will be determined separately for r < — and r > 2-. When 

< r < 



Pr ((fc — 1) relays participate) 



= (1 — Pr (a relay is in outage)) 



fc-i 



- [1-P 
= 1. 

When r > the probability that a particular relay partici- 
pates is given by 

Pr (Rj participates in second phase) 



Pr 



(£Ml + Ptol a )>(£ + e)logp) 



Pr(log(l + p| fli n> (l + -e)logp 



= Pr(b| 2 >P' 

= 

= 0. 



Therefore, when r > -2-, 

' TCI ' 

Pr (fc relays participate in second phase) = , 2 < k < n. 

Hence, when r > — , all the relays are in outage with 
probability 1. 

Consolidating the above facts, we get 

f p-("- fc )( 1 - 2 F ; ) + , < r < 
Pr(E k ) = { o , r > 2< k < n ■ (34) 

I 1 , r > k = 1 

V TCI 5 

3) Outage Probability of the NSDF Protocol: The proba- 
bility of outage of the NSDF protocol can be calculated as 
follows: 



Pout(R) 



k=l 



J2[ l_\ ) Pr(-Efe) Pi(H k in outage|f7 fc ) (35) 



It follows from ( 13~4l that for r > — no relay will participate 
in the second phase and the channel will be as shown in d32| ). 
Hence, from d33l , we can see that, 



d(r) = (1 -r)+ , — <r < 1. 

m 



(36) 



For < r < by substituting d3~TT l and (f34t in d35l ). we get 



min ^ (n - 1) 1 + (1 - r). 

2 < k < n \ \ p 



d(r) 



(n-fc) I 1 "— J 

By solving the optimization problem, we get the statement in 
Theorem [5] Please see Appendix [III] for the details. ■ 



B. DMT Optimal Codes for NSDF Protocol: 

In this subsection, we construct a DMT optimal code for the 
NSDF protocol. Once again, we shall use CDA's to construct 
ST codes and derive a code for the NSDF protocol from the 
set of matrices comprising the ST code. In Theorem [5] we 
have allowed p and q to vary with the multiplexing gain r. 
Now, for a fixed K = |, where p and q are relatively prime, 
we outline the construction of a ST code when (n — 1) relays 
are employed in a cooperative network^. 

Let t = p + nq. Consider the CDA D(L/F, a, 7), as de- 
scribed in Section ITl-D. 21 with the maximum subfield L being 
a degree t cyclic Galois extension of the field F. Consider the 
space time code X comprising of matrices corresponding to 
the left-regular representation, as in (TTTlR, of all the elements 
in the CDA D. Let Z denote the normalized code 

z = {ex I x e x} 

where 9 is chosen to ensure that 

\\9X\\ 2 F < tp, for all x ex. 

The transmitted code matrix, denoted by Z, will be of the 
form 



£ 7^t-i) 7^ 2 

£l (T(to) 

£ 2 a(h) a 2 {£ ) 



H- 

ia 2 (£t- 
,2, 



a(£ t - 2 ) <7 2 {i t -z) 



7^(4) 



7<r 



Let P out (rlogp) := p 



-d(r) 



where £i are as defined in d20l ). In the first phase the source 
transmits the first p rows of Z and in the second phase the 
source transmits the next q rows of Z. The j th relay Rj, 2 < 
j < n, if not in outage, will decode the p rows transmitted 
by the source and will transmit q rows of Z numbering from 
p + ( j — l)q +1 to p + jq. Total delay of the DMT optimal 
code will be (p + q)(p + nq) time slots. 

Theorem 6: The ST code Z described above achieves the 
DMT of the fixed-NSDF protocol. 

4 We considering only rational values of k here although, while computing 
the DMT, k was allowed to take on irrational values too. 

5 The parameter n in )17t and in this section refer to two distinct entities. 
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Proof: We will prove that the given code is DMT optimal 
irrespective of the number of participating relays. 

From [1] and [3], it follows that the given code is approx- 
imately universal and hence, DMT optimal for any channel. 
Also, from [1, Theorem 12], we know that the code Z remains 
approximately universal even if we remove a particular set of 
rows from all the matrices in the code. 

Now, during the first phase, the source transmits the first 
p rows of Z and this code will be approximately universal 
for the channel seen by the j th relay. So, when Rj is not 
in outage, it will be able to decode the source's transmission 
and participate in the second phase. Consider the case when 
k — 1 relays, where 2 < k < n, are participating in the second 
phase. The channel model will be as shown in (f29b . with Hk 
being a m x (p + kq) matrix. Since each node (either source or 
relay) transmits only q rows of the code matrix Z, the code for 
the channel in ( 1291 corresponds to deletion of (n — k)q rows 
from Z. Therefore, the row-deleted code is DMT optimal for 
the channel in d29l ) for all fc = 2, 3, . . . , n. When no relays 
participate in the second phase, the channel model is as shown 
in d32"l ). Once again, the row-deleted code is DMT optimal for 
this channel. Hence, the given code is DMT optimal for the 
NSDF protocol. ■ 

By constructing a DMT optimal code for each value of the 
ratio k = 2, we can construct DMT optimal codes for the 
variable-NSDF protocol. 

Remarks: We mention the salient features of the results in 
this section below. 

• In the class of static AF and DF protocols, the variable- 
NSDF protocol has the best DMT for any number of 
relays, except for the case of two relays. 

• The DMT of the variable-NSDF protocol for the case of 
two relays is better than the tradeoff of the SAF protocol 
[10] for r > 0.6 (see Fig. [3]). 

• For k in the range 1 < k < - J1 T , the fixed-NSDF 
protocol has a better DMT than that of the NAF protocol 
for any number of relays. 

• For k = 1, the fixed-NSDF protocol and the NAF 
protocol have the same DMT. 

. When p = q = 1, the DMT of the fixed-NSDF protocol 
coincides with that of the NAF protocol. However, the 
DMT optimal code for the fixed-NSDF protocol has a 
delay 2(n + 1), where (n—l) is the total number of 
relays, which is considerably shorter than the delay for 
the DMT optimal codes for the NAF protocol constructed 
in [9] for n > 3. The codes in [9] have delay 4(n — 1). 

• Surprisingly, for the case of one relay the optimal ratio 
k 2 turns out to be the Golden Number, k 2 



1+V5 | 
2 • 



IV. Orthogonal Selection Decode and Forward 
Protocol 

In this section, we consider the orthogonal selection decode 
and forward (OSDF) protocol. The OSDF protocol is the same 
as the NSDF protocol, except that the source remains silent in 
the second phase. 



We state the DMT of the variable-OSDF and fixed-OSDF 
protocol, but we omit the proof since the DMT can be obtained 
along similar lines to the derivation of the DMT for the 
NSDF protocol. Also, an approximately universal CDA code 
of dimension (p+(n— l)q) x (p+(n— l)q), where (n— 1) is the 
number of relays, will be DMT optimal for the OSDF protocol. 
The transmission of various rows of the code matrices by the 
source and relay will be on similar lines to that mentioned in 
Section IIIFBI for the NSDF protocol. 

A. DMT of OSDF Protocol 

Theorem 7: The DMT of the variable-OSDF protocol is 
given by 

f n fl- ( "- 1)(B " +1) < 

d(r) = 



n(l-r) 
(n-l)r+l ' 



< r < 

ft 71+1 

< r < 1 



where k„ = — ZL r. 

As in the variable-NSDF protocol, we have allowed p and 
q to vary with the multiplexing gain r. We choose the value 
of (p, q) which maximizes the diversity, for a given r, for the 
case p > q. Suppose p = nq, then the optimal value of n is 
given by 

< r < -i-r 



^n ? 

l+(n-l)r 1 
(n-l)(l-r) ' K„+l 



< r < 1 



For a fixed choice n = 2 , 
protocol is given by: 

if 1 < k < K n , 



d(r) = 

else if k > K n , 
n ( 1 



the DMT of the fixed-OSDF 




< r < 
("-i)p 



(n-l)p 
nm—p 

< r < 1 



d(r) 



m(n— 1) 
nq 



Mi - ^ 

(l-r) 



< r < ± 

— — m 

JL < r < rip-m 

m — — m{n— 1) 

np-m < < (n-l)p 

m(n— 1) — — nm—p 

(nz^P <r<i 

nm—p — — 



V. Non Orthogonal Amplify and Forward 
Protocol 

In this section, we construct a code that is DMT optimal 
for the non-orthogonal amplify-and-forward (NAF) protocol. 
For the sake of completeness, we will reproduce here the 
description of the NAF protocol. The DMT of the protocol 
was first computed in [5]. 

In this protocol, the source S transmits at each time instant, 
and the relays take turns in transmitting an amplified version of 
a previously received signal. If the number of relays is (to— 1), 
the set of equations describing a 2 (to— 1) -length frame are (see 
[5]) 



yt = 9ix t + w t , 

Vt = bihi(giX t -\ + Vi) + gix t + w t 



t odd 
= - + 1 
t even 



(37) 
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where 6j is the amplification factor at relay Ri, x t is the Let f*A^-x~ y. = (yn y# 
signal transmitted at time instant t, y t is the signal received 
at the destination at time t and Wt is the noise added at the 
destination at time t. Vi is the noise added at the relay Ri in 
this frame. 

Theorem 8: [5, Theorem 4] The DMT of the NAF protocol 
is given by 



x^BB^x = 



y ip ). Hence, 

2 

P 

E 



E^' 

3=2 



d(r) = (1 - r)+ + (n - 1)(1 - 2r)+. 

A. Explicit DMT optimal codes for NAF protocol 

We will present an explicit construction, based on CDA, 
which achieves the DMT of the NAF protocol. 

Consider a 2(n - 1) x 2(n - 1) DMT optimal CDA ST 
code. In accordance with the NAF protocol in [5], let the 
source continuously transmit the vector [x\, X2, • • ■ , £4(n-i) 2 ]> 
coming from a row by row vectorization of the code. Each 
intermediate relay Ri, i = 2, • • • , n, forwards at time t = 
4(n — l)(i — 2) + 2(n — 1) + k what it received at time 
t = 4(n - l)(i - 2) + k where k = 1, 2, ■ ■ ■ , 2(n - 1) . 

Theorem 9: The above scheme achieves the DMT of the 
NAF protocol. 

Proof: For the single relay case, we use the equivalent 
representation of the channel for the NAF protocol in matrix 
form: 



+ 

/l2«2,3 + w ± 



yi 


V3 




9i " 




Xi x 3 


2/2 


Vi _ 




32^2 gi 




X 2 XA 



^■2^2.1 + W2 



On vectorizing the above channel, we can see that the noise 
vector is white in the scale of interest. Hence, the DMT of 
the above channel is met by the approximately universal 2x2 
CDA code. 

Proceeding as in the single relay case, we can show that the 
DMT of the NAF protocol is achieved by the corresponding 
approximately universal 2(n— 1) x 2(n — 1) CDA code. ■ 

The above code for the NAF protocol was first presented 
in [17]. Around the same time, in [9], the authors constructed 
DMT optimal codes for the NAF protocol which have shorter 
delays than the codes presented here. 



Appendix I 
Proof of Lemma[T| 

Proof: Let fj = g 3 hj. Then B = J2j=2 fj^j- We have, 
x^BB^x = WB^xW 2 



J=2 



By applying the Cauchy-Schwarz inequality, we get 

p / n \ / n 

jb&x < E E i^i 2 E i; 

i=l \j=2 J \j=2 

p n 



(«-i)EEfei 

i=l ]=2 
n 

("■-i)Ell^!l 2 



j=2 



(--!)Eh/; a ] 



^n 2 



J=2 



{n-l)g?Y,\f3?AiA\ 



3=2 

With this, we get the bound stated in the Lemma. ■ 
Appendix II 

Solution of the Optimization problem for the OAF 
Protocol 

We need to solve 
d(r) < inf u + (n- l)v . (38) 

(p—q)u+qmin{u,v} > p—rnr 

If p = m, we get the non-cooperative case with the DMT 
given by 

d{r) = (l-r) + . 

Now, we separately consider different ranges of r. The opti- 
mization procedure is summarized in the flowchart in Fig. [5] 

Case A r > 

m 

Since r > — , p — mr is negative. Hence, the infimum is 
obtained by choosing u = v = and 

d(r) < 0. 

But, by definition, d(r) cannot be negative. Therefore, 

p 

d(r) = , for r > — . 

m 

Case B < r < 

— — m 

As shown in Fig|5] this case further breaks up into two cases. 
Case B.I min{w, v} = u 
Consider the inequality 



pu> p — mr, 



which leads to 



u > 1 



mr 
P 
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d(r) < 'mi(p-q) u +minu,v>p-mrU+ (n - l)v 




Fig. 5. Flowchart describing the computation of the DMT of the OAF protocol. 



Therefore, by substituting in ( I38t . we get 



P 



d(r) < n 1 - 



Case B.II min{u, v} = v 

One solution to the optimization problem in d38l l is 

mr 



u = v = 1 — 



P 



We can choose 



mr 

u = 1 \-6 

P 

mr f p-q 
v = 1 o 

p \ q 

where 5 is a small positive number. In this case, 

u + (n — l)v 

mr . , . { mr r p — q 
= 1 h<5+ (n-1 1 5- — - 

p \ p q 

(n - 1) 

q 



( mr \ 






V-T) 





(39) 



Now, we need to consider two separate cases depending on 
the ratio 

Case B.II.a < 
Suppose 

1 > 0, 



we get 



V < « 



to 2n — 1 



In this case, we can substitute S = in (1391 to obtain the 
infimum. Hence, 



d(r) < n 1 - 



Case B.II.b -2- > ^-r 

m — 2n — 1 

In this case, we choose S as large as possible under the 
constraints u < 1 and w > 0. Suppose u < 1, 

mr 

1 +6 < 1 



(40) 



Suppose f > 0, we have 



1 _rnr_ s ,P_q- q 



« < I"- • (41) 

p-q \ p J 

Hence, J is chosen to meet one of the upper bounds in (l40t 
or ( HTT i depending on r. By equating the two upper bounds, 
we get 

r 



m 



So, we need to consider two different ranges of r. 
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Case B.II.b.i < r < -2- 

— — rn 

We choose 



Then. 



d(r) < nil 



< n 1^1 
Case B.II.b.ii ■£ < r < 2- 

m 

(I 



5 = 


mr 
P 


mr 


> 


P 




mr 


1 - (n - 1 


P 


(n- 


- l)mr\ 



p- g 



rig 



We choose 



P-9 



mr 
P 



Then. 



d(r) < nil 



mr 
P 



9 



< 



p-q 
P 



1 



mr 
P 



1- (n-1) 



p-g 



1 



p- g 

Summarizing the various cases 

We make some comments with respect to the flowchart in 
Fig. [5] 

• From Case A, it follows that 

d(r) = , r > — . 

m 

. Cases B.I and B.II.a tell us that for < tt ^, 

m — In— 1 ' 

, . / ?7ir\ V 

d(r) < nil , < r < — . 

\ p J m 

This upper bound improves for higher values of 
Hence, we must choose the maximum possible value of 
-2- in the permissible range which is 2 ™_ 1 ■ 

• From Cases B.I and B.II.b, we can see that for > 



2n-l ' 



d(r) < { V J ' - - m _ 

p— q \ p J 5 m — — m 

We can see that the best tradeoff occurs for the choice 
-2- = „ " i and this tradeoff coincides with the best 

m 2n — 1 

tradeoff obtained from cases B.I and B.II.a. 
• We always have the choice of avoiding relay cooperation, 
and thus of achieving d(r) = (1 — r) for any value of 
r. 

With all the above observations, we get the statement of 
Theorem [2] 



Appendix III 
Solution of the Optimization problem for the 
NSDF Protocol 

Let 

d a (r) 



(n-l)(l-^) +(l-r)+ 



in— l)m+« 

n 1 — — —r 

np 

l-r , 



and 



d b (r) 



min fn — fc) I 1 

2<fc<n \ 



< r < 



P 
m 



< r < £- 

— — rn 

*-<r<l 



+ d k (r) 



(42) 



where dk(r) is given in d3TT ). Therefore, 



d(r) =min{ d a (r),d b {r) } , < r < (43) 

TO 



Also, from y6), we know that the DMT of the fixed-NSDF 
protocol is (1 — r) + for r > — . Since the above expression 
for the DMT also evaluates to (1 — r) + for r > for the 
sake of simplicity, we can assume that the expression in d43l ) 
is the actual DMT for all r. 

Substituting the value of dk(r) in d42l) . we can see that the 
minimum occurs when k = n. Therefore, for the case p > g, 
we have 



d b (r) 



d n (r) 



(n— l)m 

n 1 — - 1 — r 

nq 

Kl-r), 



< r < 

— — rn 

< r < 1 



Depending on the choice of the ratio n = £ and the 
multiplexing gain r, either d a (r) or dh(r) will determine the 
actual DMT d(r). It can be shown that there is a critical value 
of k, which we shall denote by k„, below which d(r) = d a {r), 
i.e., 

d(r) = d a (r) , 1 < n < n n . 

In order to compute k„, we compare the curves corresponding 
to d a (r) and db(r). We can see that at k = K n , 

(n — l)m+p (n — l)m 



np 

pnm — pm 



nq 

qnm — qm + pg 



By substituting K n = | in the above equation, we get 



(n —!)«;„ — K n — (n — 1) = 0. 



Hence, 



1 + y/1 + 4(n - l) 2 



2(n-l) 

Therefore, when k < k„, the DMT is as given in d27b . 
When k > K n , both d a (r) and c?b(r) determine the DMT for 
different ranges of r. The tradeoff is as given in d28b . We can 
check that the point of intersection of d a (r) and db{r), when 
k> k u , is given by 

np — to 



(n — 2)m + p 



(44) 
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and 



d(r) 



m 
P 



1 



np 



(n — 2)m + p 



(45) 



Hence, we have computed the DMT for the fixed-NSDF 
protocol. 

Now, to get the best possible DMT in case of the variable- 
NSDF protocol, we vary k with the multiplexing gain r, i.e., 
we choose the value of k which maximizes the diversity at any 
given multiplexing gain. By comparing the DMTs in d2Tb and 
d28l , it is clear that we must choose n = n n for r < 



l 



For r > - we need to track the point of intersection of 
d a (r) and c?&(r). This point would correspond to maximum 
diversity at a certain r. By substituting k = ^ in (l44t and 
(|45t . and by eliminating k from both the equations, we get 



d(r) 



(n — 7')(1 — r) 
(n - 2)r + 1 



1 



< r < 1. 



With all the results mentioned above, we get the statement 
of Theorem [5] 
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